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Abstract 

The crossover behavior of the semi-infinite three dimensional Ising model 
is investigated by means of Pade approximant analysis of cluster variation 
method results. We give estimates for ordinary critical as well as for multi- 
critical exponents, which are in very good agreement with extensive Monte 
Carlo simulations. 
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The present paper reports on a preliminary investigation of the multicritical behavior of 
the semi-infinite three dimensional Ising model by means of a recently proposed technique |I[ 
which combines the off-critical great accuracy of the cluster variation method (CVM) @-|4j] 
with the analysis tools for critical behavior based on the well-known Pade approximants || . 

The semi-infinite Ising model (see Ref. || for a review) is described by the (reduced) 
hamiltonian 

pH = - J s s i s i ~ J Yl s k s h (!) 

{ij)cdA (kl)<£dA 

where Si = ±1 is an Ising spin at the site i of the lattice A, dA is the free surface of the lattice, 
J s , J > are the surface and bulk nearest neighbor interaction strengths, respectively and 
(3 = 1/ksT, with ks Boltzmann's constant and T absolute temperature. In the following, 
we will take A to be a simple cubic lattice and dA a (100) free surface. 

The phase diagram of this model is well-known. For A = J s / J — 1 < A c we have the 
so-called ordinary transition at the bulk critical temperature, with simultaneous (but with 
different critical exponents) disordering of bulk and surface. For A > A c we have first the 
extraordinary transition at the bulk critical temperature, at which only the bulk disorders, 
and then the surface transition, at a higher temperature, at which the surface disorders. For 
A = A c we have a multicritical point at the bulk critical temperature, at which the special 
transition takes place. The surface transition is in the two dimensional Ising universality 
class and, in the vicinity of the special point, its critical temperature has the power-law 
behavior 

T cs - T cb ~ (A - Ac)* (2) 

where is the crossover exponent associated to the multicritical point. 

Our purpose is to show that reliable estimates for the location of the multicritical point 
and the crossover exponent, as well as other critical exponents, can be obtained by means of 
a new technique JTJ which relies on a Pade approximant analysis of CVM results. The basic 
idea is to find a region, far enough from the critical point, where the quantity of interest can 
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be obtained with a given accuracy (say, 10~ 5 ) from the CVM, and then to construct [L, M] 
Pade approximants for functions like that in Eq. |ona suitably chosen set of L + M + 1 
points in the region of accuracy (see Ref. |1[ for more details). 

In what follows, we shall use the 4 x N approximation which we have introduced in Ref. 
0. Here the semi-infinite lattice is approximated by a film of N layers, with the topmost 
layer representing the free surface, and the bottom layer constrained to the bulk, which in 
turn is studied in the cube approximation. In this system the maximal clusters for the CVM 
are chosen as those clusters with 4iV sites, formed by a column of N — 1 elementary cubes. 
The approximation is formally similar to the square (Kramers- Wannier) approximation for 
the square lattice and in fact reduces to it for vanishing bulk interaction. 

By means of the 4 x iV CVM approximation with N = 4 (but we had to use also N = 5 
for A < and low J), we determined the behavior of the surface layer magnetization ra\ 
as a function of the temperature for J > J min = 0.30 and several values of A, in the range 
—0.50 < A < 0.60. In this region the maximum error on mi is estimated |1],[7| to be about 
10~ 5 or less. 

We then determined, following the procedure proposed in Ref. [jj, Pade approximants 
for the function 

ft(z,A) = (z e -z)±lam 1 (z,A), (3) 

where z = e~ J and z c ~ 0.801 [0,|J. We obtained very stable Pade tables for j3f s (A) = 
(3*(z c , A), and our estimates are plotted in Fig. 0, which shows very good agreement with 
available Monte Carlo data 0. As already mentioned in Ref. [|J the exponent for the 
ordinary transition is estimated as (3\ ~ 0.78, again in very good agreement with Monte 
Carlo estimates P,p^|. 

For A > 0, the effective exponent is appreciably different from its limiting value fli, 
indicating that the crossover region is quite large. This allows us to study the multicritical 
behavior of our model by analyzing the behavior of the critical amplitude -Bi(A) defined by 

m 1 (z,A)~B 1 (A)(z c -zf\ z^z~. (4) 

3 



In fact, since near the multicritical point one has 

m 1 (z, A) ~ (z c - z)^rhi(x), x = , (5) 

the critical amplitude -Bi(A) must diverge as (A c — A) - ^ 1- ^/^. By fitting this behavior 
we found A c ~ 0.515 and (ft ~ PD/<P ^ 1-39, which, using /3f = /3f (A c ) ~ 0.19, yields 

~ 0.42. The calculated amplitudes, together with the theoretical curve, are shown in Fig. 
^|, where it can be seen that the fit is remarkably accurate. 

By comparing these results with those from extensive Monte Carlo simulations in Tab. 

1 we see that our values for A c and fif 1 are again in very good agreement with previous 
results, while for the crossover exponent <f), where the situation is less clear, we obtained a 
value which turns out to be closer to Ruge et al. result than to Binder and Landau one. 
As a general remark, our technique seems to be fairly more accurate than the second order 
e-expansion. It is also quite simple and does not require very large CPU times (~ 25 hours 
on a DEC Alpha for this work). 
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TABLES 



TABLE I. Comparison with Monte Carlo simulations and e-expansion. 



Present work 0.515 0.19 0.42 

Landau and Binder [9] 0.52 0.18 0.461 

Ruge et al. [10] 0.5004 0.237 0.59 

e-expansion [11] 0.245 0.68 
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FIGURES 

FIG. 1. Effective exponent for the surface layer magnetization. Present work (•, solid line is a 
guide to the eye) and Monte Carlo (x, from Ref. [9]). 

FIG. 2. Critical amplitude of the surface layer magnetization. The solid line is our fit. 
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